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■^ ■ A set of equations representing the W/Z mass ratio and various observ- 

ables of Z decays in terms of a = a{mz), G^, mz, rrit, mn, a^ = «s("^z), 
rrib and rrir (all other fermion masses being neglected) are compared with 
the latest data of the four LEP detectors, which at the level of one stan- 
dard deviation coincide with their Born values. Our global fit gives: rrit = 
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Oh! 

^_ ■ 16lllit22 , «s = 0.119 ± 0.006 ± 0.002, where the central values corre- 

spond to rriH = 300 GeV, the first errors are statistical and the second ones 

k> ! represent shifts of the central values corresponding to itlh = 1000 GeV(+) 

H I and 60 GeV(— ). The predicted mass of the top is smaller than in the recent 

fits by 4 GeV. The predicted values of rriw/fnz and the LEP observables, 
based on the fitted values of rrit and dg, show a weak dependence on ttih and 
differ by several predicted standard deviations from the corresponding Born 
values. The uncertainties of the predicted values and their deviations from 
the corresponding Born values determine the experimental accuracy required 
to observe electroweak radiative corrections. 
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Precision measurements of 5 million events of Z decays observed at LEP 
has allowed a test of the Standard Model of unprecedented accuracy. It 
turned out |I[] that all these data, for each observable separately, can be 
described within \a by a simple Born approximation, using as a basis the 
value of the running electromagnetic a at g^ = ml: 

a = (x{m\) = 1/128.87(10) , (1) 

the Fermi coupling constant: 

Gf, = 1.16639(2) X 10-^GeV"2 , (2) 

and the mass of the Z boson: 

m^ = 91.187(7) GeV . (3) 

The non-observation of genuinely electroweak radiative corrections by 
the high precision LEP experiments allows us within the framework of the 
Standard Model, to obtain rather accurate predictions of the mass of the 
top quark and of the gluon coupling constant at g^ = -ml, a^ = Q;s(m|). In 
this paper we present the results of a global fit of all LEP data and uiw/ttiz 
ratio within the Standard Model one-loop approximation based on three most 
accurately measured observables: G^, mz and a. Our results for rrit and as 
are close to but slightly different from those obtained by other authors, who 
based their global fits on other parametrizations in which the running of a 
was not separated from genuinely electroweak corrections. In particular, our 
central value of rrit is 2-4 GeV lower than in the previous fits. 

The predicted theoretical values of mw /mz and of various LEP ob- 
servables (including their uncertainties) calculated with the fitted values of 
uit ± 5mt and as ± 5as are compiled in Table 1. It can be seen that the 
predicted values differ by several standard deviations from their a Born val- 
ues. Table 1 shows what accuracy should be reached by experiments in order 
to observe the electroweak radiative corrections and to test in this way the 
Standard Model. 

Four relations describing four "gluon-free" observables miy/m^, g^, qa 
and gv I Qa in terms of a = a{rn?z)^ G^ and m^, and of the masses of the top 
quark mt and Higgs boson itih have been derived in ref. in electroweak 
one-loop approximation: 

3c 

niw/niz = c + « 3^^^2(c2 _ ^2) ^"^(^' ^) ' (4) 



^^ = 4-^64^^^^^''^^ ' (^) 



9v/9A = l-As' + a^-^^—^^Vnit,h) , (7) 

where t = {rrit/mz)'^, h = {mH/mzY, while c = cos 6* and s = sin 6* were 
defined by 

From Eqs. d), (|), (D, (D it follows that 

s^ = 0.23118(33) . (9) 

Explicit expressions of Kn, Vj,, Va, Vr in terms of t and h and parameters 
s and c are given by eqs. (40), (68), (78) and (103) of ref. [Q]. 
Each function V^ is a sum of four terms: 

Vit, h)=t + Tiit) + H,{h) + Ci . (10) 

Explicit expressions for Tj(t) in terms oft, c, s are given in ref. 0: eq.(42) 
for Tm] eq. (84) for T^; eq. (104) for Tr. 

As for the Hi{h) they have been presented in by explicit functions of 
h, c, s plus certain, rather cumbersome, combinations of c and s that were 
not written out but were evaluated in g] numerically for s^ = 0.2315(3) and 
rriz = 91.175 GeV. For the updated values of the above combinations and 
constants Ci in eq. (|Tnp, see Appendix B of ref. 0]. 

In ref. 0] the contribution of virtual gluons in the quark loops have 
been taken into account; for the sum of light and heavy quark loops up to 
terms O(^) [misprints of the preprint [^ has been corrected in the text to 
be published in Yadernaya Fizika 0]]: 

, (11) 
(12) 

• (13) 
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The expressions ([TT| ) - ([l3|) are valid for t > 1. For rrit < rriz we either put 
SVm = , 6Va = , 6Vr = 0, or used interpolation to the massless limit in 
which 

SV^ = —4{c^ - s^) Inc^ = -^ ■ 0.57 (14) 

TT vr 



5Va = ^4 (c' -s' + ^/) = +^ ■ 2.63 (15) 

SVn = 0. (16) 

As the corrections are anyway small at t < 1 they practically do not 
change the fitted values of rrit and as- In the global fit condition, we also 
used t > 1, which is definitely valid in view of the recent results of CDF 
and DO (m^ > 108 GeV at 95% CL). This latter procedure does not change 
the fitted values of rrit and a^. 

In ref. the same approach based on the a, Gf^, rriz parametrization 
and on the functions Vi was used to calculate the hadronic (qq) decays of Z. 
In the gluonless approximation the differences of V^'s for quarks and leptons 
are given by eqs. (7)-(10) of ref. |^. In this approximation the largest 
two-loop terms should be considered, which are proportional to mf: one of 
them is in the ti contribution to the Z boson self energy, the second - in the 
ti contribution to Z —>■ bb vertex. These terms have been calculated in ref. 
[0]. They become noticeable for large values oirriH/rnt. To take into account 
the first of them one has to multiply the term t in eq. (|T0|) of the present 
paper by factors in brackets derived from equations (16a) and (17a) of ref. 
[0]. To take into account the second term one has to multiply the term t in 
the eq. (18) of ref. by factors in brackets from equation (16b) and (17b) 
of ref. |7[. Note that Eqs. (16) refer to nin = 60 GeV, while Eqs. (17) - to 
tuh = 300 and 1000 GeV. 

In order to take into account the virtual gluons in the Z ^ bb vertex one 
has also to multiply the term t in eq. (18) of ref. by a factor (1 — 2.29as/7r) 
calculated in refs. [§1,0. AH other gluonic corrections, up to (ag/vr)^, known 
in the literature, are included in the eqs. of ref. 0. In particular the running 
mass of the 6-quark was also included in ref. M . We present the results of our 



fit for mb{mz) = 3.1 GeV according to refs. |]T0[, |[Tl|. If we change m},{mz) 
from 2.4 to 3.4 GeV, the predicted value of Rb at rnn = 300 GeV in Table 1 
decreases from 0.2164 to 0.2158, the fitted central value of nit decreases by 
0.5 GeV, while that of as increases by 0.002. The latter will change the Born 
values of F/j, Tz, Ch, Ri- AH masses of fermions lighter than b give very small 
contributions. Thus inclusion of m,- decreases F^ by 0.19 MeV. As for the 
running mass of the charmed quark {mc{mz) < 1 GeV), we neglect it. 

It should be emphasized that not all two-loop corrections of the order of 
aas have been calculated in the literature. In particular the vertex triangle 
electroweak diagrams with a gluon connecting a quark line in the triangle 
with an external quark line are unknown. Such two-loop corrections may 
substantially change the values of rrit given in Fig. 3, but its overall fitted 
value is expected to be changed by approximately 1 GeV. This can be seen 
by using two different expressions for F;,: one in which the specific correc- 
tions due to the ttW triangle and those due to external gluons are treated 



separately for the vector and axial channels (ref. 0, eq. (27) [|), the other 
in which they enter as a single correction (ref. §], eq. (17)). 

The equations of refs. [Hl-Pl, M - W3 described above and the latest 



experimental LEP data |T^, [^ (see Fig. 3 and Table 1) are used in this 
paper to fit the values of rrit and a^. Assuming mH= 60, 300 and 1000 GeV, 
the results of the fit are: 

mt = 162tl',tll (17) 

a, = 0.119 ±0.006 ±0.002 , (18) 

xVd.oI. = 3.5/8 , (19) 

where the central values correspond to rriH = 300 GeV, the first error is 
experimental, the second one corresponds to the variation of mn'- sign + 
corresponds to mn = 1000 GeV, sign — to tjih = 60 GeV. (Note that 
according to ref. |T^, for mn = 60 GeV and rrit from 110 to 190 GeV the 



vacuum is unstable but with a lifetime > 10^° yr.; for rrit > 190 GeV the 
vacuum is dangerously unstable. For mn > 300 GeV vacuum is stable for 
any reasonable value of m^) 

Independent constraints on rrit are given by the measurements of the 
mw/tTiz mass ratio on pp colliders by UA2 |T5| and CDF [jT6[ experiments. 



The Particle Data Group fit of the mw/rnz ratio from these experiments |]T7 
is: 

mw/mz = 0.8798 ± 0.0028. (20) 

The combined fit of LEP and pp collider data gives: 

mt = IQltlltlt , (21) 

«, = 0.119 ±0.006 ±0.002 , (22) 

xVd-o.f. = 3.5/9 . (23) 

The allowed region of rrit and as is shown in Fig. 1. The value of the 
Higgs mass is fixed here to mn = 300 GeV. The lines represent the s-standard 
deviation ellipsoids (s=l, 2, 3, 4, 5), corresponding to constant values of x^ 
(x^ = Xmin + ■5^)- So in the case of validity of Gaussian errors the projections 
of the ellipsoids on the rrit and ag axis define s-standard deviation confidence 
intervals of the corresponding parameters. 

The existing data do not put real limits on the mass of the Higgs. The 
corresponding contour plots are shown in Fig. 2, where as was fixed to as = 
0.119. The minimum of x^ is at a very low value of mn, but even the two- 
standard deviation contour is not contained in the 1000 GeV range of Higgs 



mass. Similar results were obtained by many authors [0, |T3|, |T8|-|2T| by 



^Note a misprint in the sign of the third term in the third hne of eq. (27) in ref. 



using various computer codes |22[-|2^. Our fit favours sliglitly smaller {^ 4 
GeV) values of mj. 

The individual values of rrit for each of the LEP observables for ttih = 
SOOGeV and ds = 0.119, are shown in Fig. 3. The data on A^, Ri and ahad 
are insensitive to the top mass. The data on Rf, and Ri are compatible with 
low values of the rrit, excluded by the direct search on the Tevatron. 



According to ref. |BD[ there should exist substantial non-pert urbative 



parts of gluonic corrections, which appear at the tt threshold and enter the Z- 
boson propagator through dispersion relations. These non-perturbative cor- 
rections, according to ref. [|30[, may reach 25-50% of the perturbative ones. 
The authors of ref. |^ insist that these corrections are fully non-controllable 
and eventually would prevent extraction of any information about the Higgs 
from LEP precision measurements. We consider these results as artifacts of 
dispersion calculation of Feynman integrals. When real parts of the same 
integrals are calculated directly (without referring to their imaginary parts), 
it is obvious that the non-perturbative effects are absolutely negligible, be- 
cause of the asymptotic freedom of QCD and large virtuality of ti loop for 
external momenta not larger than mz '■ 2mt — mz > 200GeV ^ ^qcd- (See 



for instance ref. |32], where these arguments were applied to charmonium.) 
Thus we neglect tt threshold effects in our fit. 

However we should note that another source of uncertainty in 0{aas) cor- 
rections really exists - the virtuality of quarks in the vector boson self energy 
loops varies with varying internal and external momenta, which influence the 
value oi as [§| . 

It is instructive to use the fitted value of rrit in order to predict the theo- 
retical values of every LEP observable for three values of rriH = 60 , 300 and 
1000 GeV. The corresponding central values of mt are, according to Eq. (p!7|), 
rrit = 139 , 162 and 178 GeV, respectively. The results of this procedure are 
presented in Table 1. One can easily see that the central predicted values of 
observables depend on rriH rather weakly. One can also see that the uncer- 
tainties of the predicted values of observables are much smaller than their 
present experimental uncertainties; also, the predicted theoretical central val- 
ues differ from the corresponding Born values by several standard deviations, 
calculated from the fitted errors of rrit and as- For instance, in the case of Rb 
it is 4.5(T-6.5cr while in the case of my/ /mz it is 3cr. The latter fact was used 
recently by Sirlin |3^ to argue that in the case of rriw /mz the existence of 
radiative corrections has been proved at the 3o" level. It is obvious, however, 
that the high accuracy of the predicted theoretical value of any of the observ- 
ables in Table 1 cannot serve as evidence for the observation of electroweak 
radiative corrections. Only direct high-precision measurements could provide 
such evidence. At present the experimental uncertainties accommodate com- 
fortably both the a Born values and the central one-loop-corrected values. 



The quality of the description of the data is not very different at the 
tree level and at the one- loop level and weakly depends on rriH- Of course 
we do not think that the a Born approximation would provide a universal 
description of all LEP data when their accuracy would be much better than 
now. However at present it seems to be sufficient. 

As for the weak dependence on mn, one has to bear in mind that it partly 
follows from the correlation between mn and rrit in the fit: the heavier is the 
Higgs the heavier is the top. Therefore when the top is discovered and its 
mass is known the limits on the Higgs mass will be more stringent (see 
columns 4-7 of the table in ref. [|I|). 

Inspection of Table 1 reveals that before the top is discovered the natu- 
ral experimental strategy is to reduce the experimental uncertainties of the 
observables to the level of the predicted ones (by a factor of 4 for Rb and 
cT/i, by a factor of 2.5 for niy^ /niz)- After the mass of the top is measured 
with accuracy ±5 GeV the reduction of uncertainties for gv/ Qa by a factor 
of 4 will be especially promising. Such a program requires not only increas- 
ing the LEP statistics by a factor of 20 (10*^ Z bosons), but - what may be 
much more difficult - reaching a qualitatively new level of sophistication in 
controlling systematics. If this can be done, LEPl will test the electroweak 
corrections of the Standard Model and may possibly reveal new physics. 

We do not discuss in this paper other possible manifestations of elec- 
troweak loops, such as in z/e scattering, deep inelastic vN interaction, parity 
violation in atoms, K^ <-^ K^ , B^ *-^ B^, B^ *-* B^ transitions, radiative 
corrections to superallowed j3 decays, such as ^^O. In some of these experi- 
ments the accuracy is much worse than in Z° decays, in others - the results 
of calculations heavily depend on nonperturbative QCD effects which can be 
estimated only roughly. As for ^"^O the virtual W boson serves here only to 
cut off the logarithmically divergent large electromagnetic correction. The 
genuine electroweak corrections must be isolated from the trivial electromag- 
netic ones in order to see whether the former can be quantitatively compared 
with available experimental data on ^'^O. 

Finally, we would like to present in Table 2 our fit in the traditional 
form used by the LEP collaborations in spite of its obvious shortcomings as 
compared with Table 1: it contains only two of six independent observables 
{syy = sin^6'vy = 1 - m^/ml and s^f^ = siia^ O^^f} = |(1 - gv/gA)) and 
the way it averages the predicted and measured values of s^ is potentially 
misleading. 

The second and the fourth columns of Table 2 are taken from ref. |13|, the 



third and the fifth present similar results of our fit taken from our Table 1. 
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Table 1 



Observable 


Exp. value 


Predicted values for observables, a Born values 
and radiative corrections. 


tuh = 60 GeV 


ruH = 300 GeV 


ruH = 1000 GeV 


mt = 139(17) 
GeV 


mt = 161(16) 
GeV 


mt = 178(14) 
GeV 


as = 0.116(6) 


a, = 0.119(6) 


a, = 0.121(6) 


mw/mz 


0.8798(28) 


0.8796(11) 


0.8799(10) 


0.8799(9) 


0.8768(2) 


0.8768(2) 


0.8768(2) 


0.0028(11) 


0.0030(10) 


0.0031(10) 


9A 


-0.50093(82) 


-0.50067(37) 


-0.50092(33) 


-0.50096(30) 


-0.50000(0) 


-0.50000(0) 


-0.50000(0) 


-0.00067(32) 


-0.00092(33) 


-0.000096(32) 


9v/9a 


0.0716(28) 


0.0711(21) 


0.0704(19) 


0.0697(17) 


0.0753(12) 


0.0753(12) 


0.0753(12) 


-0.0042(19) 


-0.0049(19) 


-0.0056(18) 


Ti (GeV) 


0.08382(27) 


0.08373(15) 


0.08381(13) 


0.08382(12) 


0.08357(2) 


0.08357(2) 


0.08357(2) 


0.00017(13) 


0.00024(13) 


0.00025(13) 


r, (GeV) 


1.7403(59) 


1.7374(40) 


1.7384(39) 


1.7382(38) 


1.7394(36) 


1.7407(36) 


1.7420(36) 


-0.0020(23) 


-0.0023(24) 


-0.0038(22) 


Tz (GeV) 


2.4890(70) 


2.4888(48) 


2.4905(45) 


2.4904(43) 


2.4877(36) 


2.4890(36) 


2.4903(36) 


0.0011(33) 


0.0015(34) 


0.0001(33) 


ah (nb) 


41.560(140) 


41.463(34) 


41.466(34) 


41.468(34) 


41.462(34) 


41.450(34) 


41.438(34) 


0.000(9) 


+0.016(9) 


+0.027(9) 


Ri 


20.763(49) 


20.749(41) 


20.742(41) 


20.738(42) 


20.815(41) 


20.830(43) 


20.845(43) 


-0.066(5) 


-0.088(5) 


-0.107(5) 


Rh 


0.2200(27) 


0.2168(6) 


0.2160(5) 


0.2154(5) 


0.2197(0) 


0.2197(0) 


0.2197(0) 


-0.0029(5) 


-0.0036(5) 


-0.0042(5) 
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Table 2 



Parameter 


LEP 


LEP + collider and vN data 


ADLO |TB 


Our fit 


Ref. [U\ 


Our fit 


mt (GeV) 
as 


1 (S(S+l'^+^9 
iDD_i9_22 

0.120(6)1^ 


162lf,li 
0.119(6)1^ 


1 r.A+l&+l% 

-L 04-17-21 

0.120(6)1^ 


0.119(6)1^ 


xVd.o.f. 


3.5/8 


3.5/8 


4.4/11 


3.5/10 


■mw (GeV) 


0.2324(5)1^ 

0.2255(19)+^ 

80.25(10)1^ 


0.2324(6)1^ 

0.2258(19);^ 

80.23(10)1^ 


0.2325(5)1^ 

0.2257(17);^ 

80.24(9)1^ 


0.2324(5)1^ 

0.2258(16);^ 

80.23(8)1^ 
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Figure and Table Captions 

Fig. 1: Allowed region of nit and ag with rriH = 300 GeV. The lines 
represent the s-standard deviation ellipsoids (s=l, 2, 3, 4, 5), corresponding 
to the constant values of x^ (x^ = xLin + ^'^)- 

Fig. 2: Allowed region of rrit and rriH with a^ = 0.119. The lines repre- 
sent the s-standard deviation "ellipsoids" (s=l, 2, 3, 4, 5), corresponding to 
the constant values of x^ (x^ = Xmm + '^^)- 

Fig. 3: The fitted values of rrit from the individual observables measured 
at LEP and pp colliders, assuming rnn = 300 GeV and a^ = 0.118. The 
hatched region corresponds to rrit < ^^z, which is definitely excluded by the 
direct search for the top quark. The central values of Rb and Ri are in the 
excluded region. 

Table 1: Observables: their experimental and predicted values. The first 
column contains observables (four of them - qa , Qvl 9a , T/ , T/i ^ are sec- 
ondary ones connected by well-known phenomenological relations with the 
observables of Fig. 3). The second column contains the experimental value 
of the observables. Columns 3, 4, 5 have three lines for each observable: the 
first line contains the predicted values which are the sum of the a Born and 
one-loop contributions for mij = 60 , 300, and 1000 GeV. The second line 
contains the a Born values, while the third gives one-loop contributions. The 
values in columns 3,4,5 are calculated for the fitted values of rrit ± Srrit and 
(is ± Sag, for each value of mn respectively. Figures for uncertainties are 
underlined when the coefficient in front of 6as or ^^^t is negative. 

Table 2: Comparison of our fit with that of the LEP Collaborations and 
the LEP Electroweak Working Group (see Table 24 of ref. |13|); sin^ ^^Tr = 
^{l—gv/gA), s"^ = sin^ 9w = ^—^TT-w/'m^'z- The collider and uN data referred 
to in columns 4 and 5 were obtained in the collider experiments UA2 |]TB| and 



CDF [|ig (see Eq. (20)) and the uN experiments CDHS 0, CHARM H 
and CCFR §^. The latter been fitted (see ref. g^) by s^ = 0.2256(47) 
{mw/mz = 0.8800(27); mw = 80.24(25) GeV). The central values in the 
table correspond to mn = 300 GeV, the upper and lower shifts to mn = 1000 
GeV and 60 GeV respectively. 
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Figure 1: Allowed region of rrit and (is with itlh = 300 GeV. The lines 
represent the s-standard deviation ellipsoids (<s=l, 2, 3, 4, 5), corresponding 
to the constant values of x^ (x^ = xLin + ■5^)- 
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Figure 2: Allowed region of mt and itlh with as = 0.119. The lines represent 
the s-standard deviation "ellipsoids" (s=l, 2, 3, 4, 5), corresponding to the 
constant values of x^ ix^ = xliin + s^)- 
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L/N 



Q 



FB 



^F 
A^F 

a: 

A. 

A' 
R, 



FB 



CT 



had 



0.7/56 ± 0.0047 
0.2320 ± 0.0016 
0.075 ± 0.015 
0.099 ± 0.006 
0.130 ± 0.025 
0.139 ± 0.014 

z 0.8798 ± 0.0028 
0.2200 ± 0.0027 
0.0158 ± 0.0018 
20.763 ± 0.049 

J 41.56 ± 0.14 nb 
2.489 ± 0.007 GeV 
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Figure 3: Value for mt from different experimental measurements, ag was 
fixed to 0.119. 
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